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ABSTRACT 

Context. Transient short-period (< 100 s) oscillations have been found in the X-ray light curves of three novae during their super-soft 
source (SSS) phase and in one persistent SSS. 

Aims. We pursue an observational approach to determine possible driving mechanisms and relations to fundamental system parame¬ 
ters such as the white dwarf mass. 

Methods. We performed a systematic search for short-period oscillations in all available XMM-Newton* and Chandra** X-ray light 
curves of persistent SSS and novae during their SSS phase. To study time evolution, we divided each light curve into short time- 
segments and computed power spectra. We then constructed a dynamic power spectrum from which we identified transient periodic 
signals even when only present for a short time. We base our confidence levels on simulations of false-alarm probability for the chosen 
oversampling rate of 16, corrected for multiple testing based on the number of time segments. From all time segments of each system, 
we computed fractions of time when periodic signals were detected. 

Results. In addition to the previously known systems with short-period oscillations, RS Oph (35 s), KT Eri (35 s), V339 Del (54 s), 
and Cal 83 (67 s), we found one additional system, LMC 2009a (33 s), and also confirm the 35 s period from Chandra data of KT Eri. 
The oscillation amplitudes are of about < 15% of the respective count rates and vary without any clear dependence on the X-ray count 
rate. The fractions of the time when the respective periods were detected at 2cr significance (duty cycle) are 11.3%, 38.8%, 16.9%, 
49.2%, and 18.7% for LMC 2009a, RSOph, KTEri, V339Del, and Cal 83, respectively. The respective highest duty cycles found in 
a single observation are 38.1%, 74.5%, 61.4%, 67.8%, and 61.8%. 

Conclusions. Since fast rotation periods of the white dwarfs as origin of these transient oscillations are speculative, we concentrate 
on pulsation mechanisms. We present initial considerations predicting the oscillation period to scale linearly with the white dwarf 
radius (and thus mass), weakly with the pressure at the base, and luminosity. Estimates of the size of the white dwarf could be useful 
for determining whether these systems are more massive than typical white dwarfs, and thus whether they are growing from accretion 
over time. Signs of such mass growth may have implications for whether some of these systems are attractive as Type la supernova 
progenitors. 

Key words, novae, cataclysmic variables - stars: individual (V339 Del) - stars: individual (RS Oph) - stars: individual (KT Eri) - 
stars: individual (LMC 2009a) - stars: individual (V2491 Cyg) - stars: individual (V4743 Sgr) - stars: individual (LMC 2012) - stars: 
individual (V959 Mon) - stars: individual (V1494 Aql) - stars: individual (V5116 Sgr) - stars: individual (U Sco) - stars: individual 
(HV Cet) 


1. Introduction 

The class of super-soft sources (SSS) was empirically de¬ 
fined as X-ray sources emitting a blackbody-like spectrum 
yielding temperatures between 20-100eV (3 -7x10^ K) and 
luminosities above 10^^ erg s ' (lGreineilll996l) . The first such 


* XMM-Newton is an ESA science mission with instruments and con¬ 
tributions directly funded by ESA Member States and NASA. 

** Software provided by the Chandra X-ray Center (CXC) in the ap¬ 
plication package CIAO was used to obtain science data. 


source s were found with the Einstein Observatory bv iLong et al.l 
(Il98ll) and were later defined as a class after more were 
discovered with ROSAT (iTruempedl 1 9921: iKahabka et aljfl9^ 
iGreiner et al.l 1199 111 . It is now commonly accepted that SSS 
emission originates in binary systems containing a white dwarf 
primary that hosts nuclear b urning of material tha t is acc reted 
from a secondary star; see Ivan den Heuvel et alJ (1 19921) and 
Kababka & van den HeuveH (Il997l) . Based on this interpretation, 
the term SSB (super-soft X-ray binary) is sometimes used. In 
high spectral resolution data that can be achieved with the 
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gratings on board XMM-Newton and Chandra, details can be 
seen that clearly confirm theoretical expectations that the X-ray 
spectra are not a blackbody (as e.g. several neutron stars), 
but an atmospheric continuum spectrum with absorption lines 
and edges. Important system parameters can be derived using 
atmsophere models, which was, for example, demons trated for 
the famous persistent SSS Cal 83 bv iLanz et al.l (l2005h . 

In addition to the persistent SSS, classical and recurrent 
novae (CNe, RNe) pass through a phase during which they 
emit an SSS spectrum in X-rays that is thought to originate 
from the hottest layers closest to the surface of the white dwarf. 
The SSS phase lasts until the hydrogen content of the accreted 
material is consumed or ejected. The X-ray and UV evolution 
of novae during their outbursts has been determined using Swift 
monitoring observations. Deeper, continuous XMM-Newton 
and Chandra observations have been obtained, guided by the 
long-term evolution of UV/X-ray emission determined by the 
Swift UVOT and XRT instruments. The turn-on time of the SSS 
phase, when bright SSS emission becomes visible, depends 
on the evolution of the nova ejecta. Accurate predictions are 
currently not possible from evolutionary models, but empirical 
scaling relations have recently been derived fr om population 
studies dHenze et alj|20l4; ISchwarz et al.ll20Tl1) . Their predic¬ 
tive power is not yet well established, and Swift is still needed 
to study the X-ray evolution of individual novae and to guide 
deeper X-ray observations. 

The hrst Swift X-ray monitoring campaign of a nova was 
performed d uring the 2006 outbu rst of the recurrent nova 
RS Oph, and lOsborne et alJ (l2006l) found short-period oscil¬ 
lations of about 35 s in some but not all X-ray observations 
taken during the SSS phase. This transient signa l in RSOph 
was studied in more detail bv lOsborne et alJ d201 ih . The period 
was hrst seen on day 32.9 after outburst and last detected on 
day 58.8, consistent with the start of the SSS phase and the 
onset of the decline, res pectively. Indepen dent evidence for this 
period was reported by iNess et al.l d2007h in an XMM-Newton 
observation taken on day 54 after outburst. They conhrmed 
that the signal came from the SSS component and not from 
the shock emission that was pre s ent in the same observation 
at higher energies. iNelson et al.l d2008h also reported having 
seen this signal for a short time in an XMM-Newton observation 
taken on day 26.1. 

Surprisingly, a transient 35s pe riod was also seen i n 
SwiftfXKY data obtained for KTEri dBeardmore et alJ l201f)h . 
Most recently, a 54s tran sient period was found for the nova 
V339Del in SwiftIX m dBeardmore et aP l2013h and XMM- 
Newton observations dNess et al.l 2013bh . both taken during the 
SSS phase. Two Chandra observations taken on 2013 November 
9 (day 87.2) and 2013 December 6 (day 114) also contain the 
54s period (Nelson private communication, publication in 
preparation). 

Similar X-ray periods were found in probab ly non-burning 
white dwarfs, for instance, a 33s spin of AE Aqr dPatterson et al.l 
Il980l) . ~ 29s in WZ Sge dNucita et al.l 12014 and references 
therein), or few-second quasi-periodic oscillatio ns (QPOs) in a 
numb er of dwarf noyae, e.g., SS Cyg, U Gem, dCordoya et al.l 
Il984ll : re feiTed to as dwarf noya oscil l ations, DNOs, see for 
example, iMauchel d2002h : IWarner et al.l d2003ah and references 
therein. Since the X-ray spectra of these systems reflect dis¬ 
tinctly different production mechanisms from SSS spectra (e.g.. 


accretion disc or the boundary layer), we do not include these 
CV systems in our study. 

In addition to noy ae, a period of ~ 67 s has been seen by 
lOdendaal et al.l d2014ll in the prominent persistent SSS Cal 83, 
which they interpreted as driyen by the rotation period of the 
white dwarf. While this is close to the break-up period, they 
argued that the white dwarf may be spun up by accretion disc 
torques. A +3 s drift from the median yalue in this scenario 
would be more difficult to explain (see the discussion in 
Sect. l5.2l l. 

All indicators suggest that the short-period oscillations 
originate in the SSS component. This type of short-period 
oscillation may thus originate from the surface of the white 
dwarf while undergoing nuclear burning. It will be of interest to 
determine whether the 35-54s periods in the noyae are related 
to the 67s period in the persistent SSS Cal 83. If the origin is 
not the spin period of the white dwarf, it might be related to the 
interiors or the nuclear burning regions, and understanding these 
processes would giye a new diagnostic method to determine 
fundamental parameters such as the mass of the underlying 
white dwarf. It is therefore important to hnd as many systems 
as possible that host short-period oscillations, ideally coyering 
a wide range of system parameters, to identify those properties 
that driye them. 

We haye searched for short-period yariations in all systems 
showing bright SSS emission, focusing here on XMM-Newton 
and Chandra obseryations, while a similar project based on 
all SwiftpCRT data will be presented by Beardmore et al. Eor 
systems with short-period oscillations, we study the eyolution 
of the power and the period. 

We briefly describe the techniques and the obseryations used 
in this article in Sect. |2] We then describe the timing analysis in 
Sect. |3] and the results of it in Sect.|4] We discuss our results in 
Sect. |5] and summarise our hndings and conclusions in Sect.|6] 


2. Observations 


XMM-Newton and Chandra obserye from highly elliptical orbits 
allowing long uninterrupted obseryations for up to two days. 
The instrumentation on board proyides high resolution in time 
and energy using yarious combinations of CCD detectors and 
dispersiye gratings. Eor more details, we refe r to the correspond - 
ing papers about the XMM-Newton satellite (l.lansen et al.ll200l]) 
hosting the Europea n Photon Imaging Ca mera, EPIC, consisting 
of the pn detector (IStriider e t alJ 2() 0lli an d the Metal Oxide 
Semi-conductor, MOS, ( Turner et al. ^O l]). and the Refl ection 
Grating Spectrometers. RGS. dden Herder e tal.1 12001 ). The 
Chandra satellite is described by IWeisskopf et ah (2002h with 
the High Resolution Camera HRC. dMurray et ^ l2000l). the 
Low Energy Transmission Grating, LETG, ( Brinkman et al.l 
l2000t). and th e Adyanced CCD Imaging Spectrometer, ACIS 
dNousek et al.l[T98^ . 


A large sample of hi gh-resolution X-ray grating spectra 
of SSS was presented by iNess et al.l (l2013al) . where a list of 
obseryations is giyen in their Table 2, from which we selected 
the obseryations in this work. We used the light curyes from the 
contemporaneous XMM-Newton[EP\C detectors or zero order 
for Chandra transmission grating obseryations. In Tables [1] and 
|2]we list the target names, day after outburst for noyae (reference 
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day given in Table 1 of iNess et al.l2013all . observation identifier, 
detector from which light curves were extracted (see table 
footnotes and text below), start date, exposure time, and the bin 
size used for the period analysis. The following five columns 
give the larger bin size used to determine count rates given in 
the following four columns, minimum and maximum count 
rates (thus the range of variability) supplemented by median and 
mean count rates. These numbers give an idea of the brightness 
and the degree of variability. In the last column we list a code 
that indicates the detrending method (see Sect. [3]). 

For XMM-Newton observations, we extracted light curves 
from the EPIC detectors pn in Is time bins. V4743 Sgr, RS Oph, 
and V339Del were too bright for the EPIC cameras, which 
suffered excessive event losses from buffer overflow. Eor these 
novae, we used the light curves from either the less sensitive 
MOS cameras (V339Del) or from the Reflection Grating 
Spectrometers (RGS). The instrument with mode and filter used 
is given in the fourth column of Table [T] Some details of the 
instruments are also included in the corresponding graphical 
representations of the results. 

The raw data were processed with standard SAS tools of 
version 13.5. We started with a standard run using the general 
SAS tool xmmextractor, which delivers a full set of event 
files and science products using standard parameter settings, 
some of them optimised to the specific observation. We then 
inspected the science products and re-generated the EPIC and 
RGS light curves over the 0.1-10-keV range using evselect 
and epiclccorr or rgslccorr. We furthermore extracted all 
available light curves that were taken with the optical monitor 
in fast mode in time binnings of 0.5 seconds. 

All Chandra observations in our sample used the LETG. We 
extracted the light curves from the zero-order photons of the 
HRC or the ACIS (Cal 87) using the regions of zero order that 
are automatically created during extraction of grating spectra. 
The extraction was made with the task Ciao v. 4.6 dmextract 
in binning of 1 s (for HRC) and 2 s for ACIS, without energy 
filtering. 


3. Analysis 

Before searching for periodic oscillations, we detrended the 
light curves. We first fitted either a high-order polynomial to 
the light curve or smoothed it with a Gaussian kernel. The 
polynomial or smoothed curves were then subtracted from 
the light curve and the difference used for further analysis. In 
the last column of Table [T] we list the detrending mode as Pn 
when using an n-th order polynomial and Gn when using a 
Gaussian-smoothed light curve, with n reporting the full width 
at half maximum (in seconds). 

The period searches were performed using p o wer s pectra 
calculated with the method of iHorne & Baliunaa (Il986l) . The 
details are described in the appendix. Since past observations 
imply that we here investigate a transient signal, we have to 
consider the possibility that a signal might not be detectable 
in a full light-curve. To detect any sporadic appearance within 
an observation, we therefore computed power spectra from 
overlapping t - 1000s time segments in adjacent steps of 500 s. 
In all cases we tested Nf - 500 frequencies within a range of 
periods between 25 s and 100s. This range covers all periods so 


far reported and avoids contamination by low-frequency noise. 

In the Appendix we present simulations supporting the 
way we interpret the resulting powers from the power spectra 
as likelihood of detecting a periodic signal. We also discuss 
corrections for multiple testing when studying various time 
intervals from a singl e light-curve based on the conserv ative 
Bonferroni correction (lBonferronill936HDun ^Il9.59[l1961h . 

To illustrate the evolution of a periodic signal, we display 
the series of power spectra with the concept of a dynamic power 
spectrum, similar to the dyna mic spectra that were used for 
X-ray spectra by iNessI (12012h . for example. A comparison is 
shown in Eig. [1] (discussed in Sect. 14.11 1. where time evolution 
maps of the X-ray spectrum (left panel) and of the power 
spectrum (right panel) are shown. The light curve is shown with 
time running down in the right part of the left panel, turned 
around by 90° clockwise, where in this case both EPIC/pn and 
EPIC/MOS light curves of RS Oph are shown in units of count 
rates. In the right plot, the modulation amplitudes are shown 
in units of counts per second for each 500 s time bin along 
the same vertical time axis. The horizontal period axis in the 
right plot is centred on the 35 s period found in the later Swift 
and XMM-Newton observations. While we always searched 
for periods over 25-100 s, we here show only the period range 
around the significant signal. We defined the colour scheme 
based on detection likelihood using red at ~ 30-80% and blue 
at > 90%. The colour scheme and the marked confidence levels 
are based on the thresholds after correcting for multiple testing 
and oversampling. The colour-encoding can be identified from 
the colour bar and the corresponding horizontal lines in the 
upper right corner. 

We first searched for significant signals within the entire 
25-lOOs period interval in each time segment. Eor the strongest 
peak, we focused on a narrow period interval +2 s to determine 
the period from a spline fit. We then determined the amplitude 
from a sine fit to the folded light curve. Rough uncertainties in 
amplitude were estimated by observing the increases in when 
varying the amplitude from the best-fit value while keeping the 
other two parameters (period and phase) fixed at their respective 
best-fit values. 

We then calculated the uncertainty of the pulse period in 
two different ways. In the first ap proach, we com puted the 
standard deviation using Eq. (4) in iLarssonl (Il996h . using the 
standard deviation of count rates in the corresponding time 
segment, cTiot, and the amplitude of the best- fit sinusoid fitted to 
the fo lded light cu r ve; se e also Eq. (14) in iHorne & Baliuri^ 
(119861) and iKova^ (Il981h . Here we neglected the uncertainties 
in the amplitude. In the second approach, we determined the 
standard deviation of the pulse periods obtained from Monte 
Carlo simulated light curves. These light curves were obtained 
by replicating the folded light curve using the same sampling 
as the original light curve and adding the appropriate Poisson 
noise. We then calculated the standard deviation of the periods 
measured from 10,000 such realizations and assigned this 
number to the uncertainty of the period. In general, both 
approaches yield comparable results. In this work we use twice 
the standard deviation as estimates of uncertainties in period. 
Assuming Gaussian statistics, this would correspond to 95% 
confidence ranges. The respective results are illustrated with the 
short horizontal black lines in the dynamic power spectra plots 
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Table 1. Journal of X-ray observations of super-soft X-ray sources 


Target 

Day“ 

ObsID* 

Detector^ 

Start time Exp. time'' 

bin size 

Min. Median 

Mean 

Max. 

D-mode* 





(UT) 

(ks) 

(s)“ (s)^ 

count rate (counts s 

‘) 


Cal 83 


0123510101 

pn/FE“ 

2000-04-23.84 

39.8 

1 

219 

2.11 

5.65 

5.00 

7.41 

G1200 

(Sect.l431 


1900 

HRC 

2001-08-15.67 

36.4 

1 

220 

0.07 

0.14 

0.14 

0.20 

P9 



0500860201 

pn/SW/m 

2007-05-13.92 

11.5 

1 

78 

4.70 

6.29 

6.30 

8.05 

G500 



0500860301 

pn/SW/m 

2007-07-06.98 

10.5 

1 

78 

4.18 

5.81 

5.78 

7.19 

G500 



0500860401 

pn/SW/m 

2007-08-21.63 

7.47 

1 

37 

2.73 

4.45 

4.52 

6.66 

G300 



0500860501 

pn/SW/m 

2007-10-05.99 

12.5 

1 

79 

2.08 

4.30 

4.40 

7.27 

G500 



0500860601 

pn/SW/m 

2007-11-24.90 

20.0 

1 

107 

3.31 

5.72 

5.53 

7.28 

G500 



0506530201 

pn/SW/tn 

2008-03-20.02 

7.27 

1 

65 

0.24 

0.66 

0.68 

1.38 

G200 



0506530501 

pn/SW/tn 

2008-04-16.52 

4.59 

1 

42 

0.86 

1.79 

1.77 

2.77 

G200 



0506530601 

pn/SW/tn 

2008-04-17.57 

10.8 

1 

42 

0.62 

1.43 

1.54 

3.44 

G200 



0506530801 

pn/SW/tn 

2008-04-19.28 

5.47 

1 

86 

0.18 

0.36 

0.41 

1.02 

G200 



0506530901 

pn/SW/tn 

2008-04-20.94 

11.2 

1 

80 

0.07 

0.48 

0.51 

1.24 

G200 



0500860901 

pn/SW/tn 

2008-04-21.09 

12.3 

1 

139 

0.10 

0.23 

0.27 

0.68 

G300 



0506531001 

pn/SW/tn 

2008-04-21.78 

8.67 

1 

79 

0.10 

0.49 

0.52 

1.48 

G200 



0506531201 

pn/SW/tn 

2008-04-23.47 

6.97 

1 

123 

-0.06 

0.24 

0.23 

0.48 

G200 



0506531301 

pn/SW/tn 

2008-04-25.34 

9.17 

1 

72 

0.12 

0.48 

0.51 

1.25 

G200 



0506531401 

pn/SW/tn 

2008-04-29.03 

13.7 

1 

214 

0.05 

0.31 

0.31 

0.91 

G500 



0506531501 

pn/SW/tn 

2008-08-12.62 

6.47 

1 

30 

5.12 

8.19 

8.09 

10.62 

G300 



0506531601 

pn/SW/tn 

2008-09-17.47 

6.37 

1 

235 

-0.06 

0.12 

0.17 

0.60 

G200 



0506531701 

pn/SW/tn 

2009-05-30.34 

45.7 

1 

105 

5.02 

7.85 

7.79 

10.04 

P9 

Cal 87 


1896 

ACIS 

2001-08-13.84 

98.8 

2 

412 

0.04 

0.08 

0.08 

0.14 

P9 



0153250101 

pn/FF/tn 

2003-04-18.87 

76.7 

1 

163 

1.30 

2.18 

2.14 

2.71 

G1200 

HVCet 

72.2 

9970 

HRC 

2008-12-18.62 

36.8 

1 

144 

0.10 

0.24 

0.24 

0.42 

G200 

KTEri 

71.3 

12097 

HRC 

2010-01-23.91 

16.6 

1 

56 

4.21 

11.67 

13.16 

27.33 

G400 

(Sect.HS 

79.3 

12100 

HRC 

2010-01-31.94 

7.22 

1 

21 

62.27 

81.29 

80.79 

98.64 

G200 


84.6 

12101 

HRC 

2010-02-06.27 

6.55 

1 

21 

25.14 

54.55 

55.48 

87.82 

G200 


158.8 

12203 

HRC 

2010-04-21.45 

7.17 

1 

21 

89.00 

107.79 

107.58 

118.88 

G200 

EMC 2009a 

90.4 

0610000301 

pn/SW/tn 

2009-05-06.43 

37.2 

1 

104 

1.30 

4.15 

4.91 

15.53 

G1200 

(Sect.lMl 

165.0 

0610000501 

pn/SW/tn 

2009-07-20.03 

57.1 

1 

100 

26.88 

47.36 

51.08 

86.23 

G1200 


196.5 

0604590301 

pn/SW/tn 

2009-08-20.59 

31.4 

1 

100 

40.47 

49.35 

51.83 

73.75 

P9 


230.0 

0604590401 

pn/SW/tn 

2009-09-23.02 

50.7 

1 

101 

15.58 

37.32 

35.65 

46.83 

P9 

LMC2012 

31.5 

14426 

HRC 

2012-04-26.91 

21.8 

1 

35 

1.50 

2.15 

2.14 

3.32 

G200 

QRAnd 


0075 

HRC 

2000-09-28.42 

53.5 

1 

284 

0.05 

0.10 

0.10 

0.17 

P9 



0047940101 

pn/SW/tn 

2001-12-31.77 

56.7 

1 

162 

1.77 

2.24 

2.24 

2.88 

G1200 

RS Oph 

26.12 

0410180201 

pn/Ti/tn 

2006-03-10.96 

10.5 

1 

30 

4.96 

17.15 

15.12 

26.79 

G200 

(Sect.lOl 

39.7 

7296 

HRC 

2006-03-24.52 

11.6 

1 

35 

14.86 

27.41 

29.53 

46.94 

G200 


54.0 

0410180301 

RGS 

2006-04-07.88 

18.8 

5 

42 

125.29 

156.63 

155.93 

169.74 

P9 


66.9 

7297 

HRC 

2006-04-20.73 

11.0 

1 

21 

62.64 

68.56 

68.67 

74.20 

G200 


“After fret (see table 1 inlNess et aL2013ai) 

^Observation Identifiers 

“On board XMM-Newton (pn, MOS, RGS) and Chandra (HRC, ACIS). For pn, mode and filter are given as FF=Full Frame, LW=Large Window, 
SW=Small Window, Ti=Timing and tn=thin filter and m=medium filter, respectively. The first XMM-Newton observation of Cal 83 consists of 
three consecutive exposures with thin, medium, and thick filter. 

“^Exposure time (ks= 10^ s) 

“Time bin size for analysis 

^Time bin size for determination of max/min count rates 
^Detrending Mode, see text. 


such as Fig. [T] centred around the respective strongest signal. 

We determined the widths of the peaks in the power spectra 
from a Gaussian fit to the peak in the power spectrum in 
frequency space to compare them to the frequency resolution 
imposed by windowing, determined by the duration of the time 
segments; Af - I IT. The corresponding width in period space 
depends on the period p, yielding jT, thus a narrower peak 
for shorter periods. If the measured width of a significant peak 
at period p is broader than p^ jT, then the period is not constant. 


Finally, we computed the fraction of the total time that a 
period is detected above a given detection threshold, similar to 
a duty cycle, where we determined these fractions for a range 
of thresholds. For each source, we computed new time series 
of power spectra that do not overlap, and in this sample, we 
summed the integration time of all time segments that contain 
a signal above a given threshold. We computed the fractions 
for each observation separately and for all time segments from 
all observations combined. The results are described in Sect. 14.71 
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Table 2. Continuation of Table [T] 


Target 

Day“ 

ObsID* 

Detector* 

Start time Exp. time^ 

bin size 

Min. 

Median 

Mean 

Max. 

D-mode^ 





(UT) 

(ks) 

(s)“ (s)f 

Count rate (counts s 

') 


RXJ0513.9-6951 

3503 

HRC 

2003-12-24.29 

49.3 

1 

271 

0.04 

0.13 

0.13 

0.25 

P9 



0151410101 

pn/Ti/tn 

2004-04-28.74 

10.7 

1 

106 

0.85 

1.84 

2.14 

4.51 

G200 



0151412101 

pn/Ti/tn 

2004-05-02.08 

15.7 

1 

59 

14.53 

20.87 

21.07 

30.08 

G200 



0151412201 

pn/Ti/tn 

2004-05-05.81 

16.7 

1 

60 

13.03 

17.78 

17.98 

25.18 

G200 



0151412301 

pn/Ti/tn 

2004-05-09.91 

14.8 

1 

59 

10.60 

18.69 

18.83 

32.10 

G200 



0151412401 

pn/Ti/tn 

2004-05-12.69 

9.48 

1 

305 

0.14 

2.68 

2.97 

9.31 

G200 



0151412501 

pn/Ti/tn 

2004-05-16.06 

12.7 

1 

462 

0.20 

2.71 

3.59 

9.70 

G200 



0151412601 

pn/Ti/tn 

2004-05-18.35 

12.7 

1 

358 

0.18 

3.97 

3.31 

7.27 

G200 



0151412701 

pn/Ti/tn 

2004-05-26.24 

16.9 

1 

358 

0.14 

2.56 

3.16 

9.29 

G200 



0151412801 

pn/Ti/tn 

2004-05-28.27 

14.7 

1 

385 

0.27 

3.34 

3.97 

9.60 

G200 



5440 

HRC 

2005-04-20.80 

26.5 

1 

249 

0.08 

0.18 

0.18 

0.29 

P9 



5441 

HRC 

2005-04-27.95 

26.5 

1 

178 

0.24 

0.48 

0.48 

0.65 

P9 



5442 

HRC 

2005-05-03.24 

27.5 

1 

179 

0.34 

0.46 

0.46 

0.63 

P9 



5443 

HRC 

2005-05-13.82 

24.3 

1 

208 

0.13 

0.33 

0.33 

0.51 

P9 



5444 

HRC 

2005-05-19.09 

26.7 

1 

193 

0.18 

0.35 

0.35 

0.58 

P9 

USco 

18.7 

12102 

HRC 

2010-02-14.49 

25.1 

1 

88 

0.19 

0.37 

0.37 

0.56 

P9 


22.9 

0650300201 

pn/SW/tn 

2010-02-19.65 

63.3 

1 

100 

13.90 

26.35 

25.90 

37.17 

P9 


34.8 

0561580301 

pn/SW/tn 

2010-03-03.61 

62.3 

1 

300 

19.67 

28.44 

28.24 

37.28 

P9 

V1494Aql 

301.5 

2308 

HRC 

2000-09-28.29 

9.94 

1 

39 

0.29 

0.67 

0.69 

1.29 

G400 


304.7 

0072 

HRC 

2000-10-01.42 

19.8 

1 

24 

0.27 

0.74 

0.95 

6.99 

G400 

V2491 Cyg 

39.9 

0552270501 

pn/Ti/m 

2008-05-20.59 

38.0 

1 

100 

53.86 

339.37 

325.90 

564.52 

G1200 


49.6 

0552270601 

pn/Ti/m 

2008-05-30.35 

30.2 

1 

108 

34.37 

59.90 

60.41 

90.64 

G1200 

V339 Del" 

87.2 

15742 

HRC 

2013-11-09.75 

50.2 

1 

100 

39.0 

46.70 

46.8 

71.4 

G200 

(Sect.l43l> 

99.0 

0728200201 MOS2/SW/m 

2013-11-21.56 

32.4 

1 

2 

63.09 

192.76 

194.92 

318.70 

P9 


99.0 

0728200201 

RGS 

2013-11-21.56 

32.5 

11 

200 

26.84 

30.28 

30.20 

32.88 

P9 


113.9 

15743 

HRC 

2013-12-06.51 

50.2 

1 

100 

23.4 

33.1 

32.7 

55.3 

G200 

V458 Vul 

460 

0555691401 

pn/FF/m 

2008-11-10.25 

9.83 

1 

130 

0.05 

0.15 

0.15 

0.30 

G500 

V4743 Sgr 

180.4 

3775 

HRC 

2003-03-19.40 

27.1 

1 

21 

0.22 

32.76 

29.83 

53.70 

G200 


196.1 

0127720501 

RGS 

2003-04-04.93 

35.3 

11 

11 

76.62 

104.00 

104.25 

138.64 

G600 


301.9 

3776 

HRC 

2003-07-18.90 

13.7 

1 

19 

30.08 

35.70 

35.73 

42.29 

G200 


371.0 

4435 

HRC 

2003-09-25.99 

14.4 

1 

28 

18.31 

28.65 

28.36 

35.28 

G200 


526.1 

5292 

HRC 

2004-02-28.06 

13.5 

1 

29 

2.37 

3.82 

3.78 

5.11 

G200 

V5116Sgr 

609.7 

0405600201 

pn/LW/m 

2007-03-05.72 

8.80 

1 

213 

12.68 

18.83 

36.98 

77.22 

GlOO 


781.8 

7462 

HRC 

2007-08-24.80 

37.6 

1 

277 

0.13 

0.21 

0.22 

0.45 

G1200 

V723 Cas 

5482 

0652070101 

pn/SW/tn 

2010-08-26.76 

50.0 

1 

318 

0.29 

0.67 

0.66 

1.04 

G1200 


6018 

0673490101 

pn/SW/tn 

2012-02-14.02 

89.7 

1 

445 

0.05 

0.24 

0.26 

0.51 

G1200 

V959 Mon 

116.0 

15596 

HRC 

2012-12-03.82 

27.1 

1 

132 

0.04 

0.14 

0.16 

0.36 

G200 


“After fref (see Table 1 inNess et al,,2013a) 

^’Observation identifiers 

“On board XMM-Newton (pn, MOS, RGS) and Chandra (HRC, ACIS). For pn, mode and filter are given as FF=full frame, LW=large window, 
SW=small window, Ti=timing, and tn=thin filter and m=medium filter. 

“^Exposure time (ks= 10^ s) 

“Time bin size for analysis 

^Time bin size for determination of max/min count rates 
^Detrending mode, see text. 

'■Reference date for V339 Del: 2013-8-14.584 


4. Results 

We performed the analysis described in Sect. [3] for the obser¬ 
vations listed in Table [1] and discuss the results for some of the 
individual systems in the following subsections. 

4.1. RSOph 

We searched for the 3 5s period found from Swift observations by 
lOsborne et al.l (l2006l) in three XMM-Newton and two Chandra 
observations of RS Oph. We first analysed the XMM-Newton ob¬ 
servation taken on day 14 (ObsID 0410180101), which is not 


listed in Table[I]b ecause it was not taken during the SSS phase; 
dNess et alJ2009t) . We divided the 23ks observation into 46 over¬ 
lapping 1000s time segments but found no significant detection. 
While a 35s signal with a false-alarm probability of 10% is found 
in one time interval, the correction for multiple testing reduces 
the detection probability to < 1%. Thus, repeating a search 46 
times in a light curve without a signal could produce such a peak 
and it may thus be a random occurrence. The XMM-Newton ob¬ 
servation taken on day 26.1 is also dominated by the early shock 
emission, but has an a dditional soft comp onent that appeared 
during the observation. iNelson et al.l (l2008ll reported that a 35s 
period was present during a short time when the soft count rate 
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Fig. 1. RS Oph on day 26.1: Spectral (left) and period (right) time- maps on the same ti me scales in vertical direction. Left: Spectral 
time-map showing the evolution of the RGS spectrum; taken from lNess et alJ (l2013ah . adding unidentified emission features with 
their wavelengths in A. The simultaneous EPIC/pn and MOS2 light curves are shown in the right panel, rotated clockwise by 90°. 
The pn has a longer initial overhead such that the exposure started ~ 1170 s later than the RGS and the MOS exposures, thus the 
black curve (pn) starts at 0.33 h, while the red curve (MOS) starts at Oh. Right: Dynamic power spectrum consisting of the power 
spectra computed from 24 overlapping 1000s time segments from the pn light curve. To map them into the 1 Iks time frame of the 
observation, the width of each row in the central panel is reduced to 500s. The horizontal orange/blue dashed lines running across 
the main panel border two adjacent 500s time intervals that correspond to the (wider) coloured bars in the right panel. They mark the 
overlapping 1000s time intervals from which the individual power spectra were computed that are shown with the same colours in 
the top panel. The colour scheme for the main panel is decoded in the top right part; blue colours were chosen to mark periods with 
>90% confidence and red for 30-80%, while yellow/light green reflect noise. Colour-coding and thresholds have been corrected for 
oversampling and multiple testing (see text). The short horizontal black line in the main panel at 35s represents the 90% uncertainty 
range in period of a significant (> 2cr) signal. In the top panel, in addition to the two power spectra from partial light-curves, the 
power spectrum from the full light-curve is shown in black; the same oversampling factor of 16.67 has been applied and contains no 
significant signal. The widths of the peaks in the respective power spectra are consistent with the intrinsic period resolution imposed 
by windowing. At the right, we show the evolution of the modulation amplitude, most of them are upper limits. 


rapidly increased by a factor 3, and in Fig.[T] we give a graphi¬ 
cal illustration of the spectral (left) and period (right) evolution 
on the same vertical time axis. Numerous s hort-lived unidenti - 
fied soft spectral emission features, noted bv iNess et alJ (l2013ah . 
are included with the values of wavelengths of unidentified fea¬ 
tures in the top left panel. So far, no convincing identifications 
were found after inv estigating possibilities of extrem e line shifts 
( Nelson et al.ll2008l) or exotic elemental abundances (lOrio et alJ 

1201 nh 

In two time segments around the same time, a periodic 
signal appeared at 35 s above 90% detection probability, which 
can be seen in the right panel. The black line in the top panel is 
the power spectrum from the full light-curve and contains no 
significant signal at 35s. While this might seem like a marginal 
result, we emphasise that the confidence levels and colour 
scheme have been conservatively corrected for the effects of 
oversampling and multiple testing with 24 trials. Moreover, 
a similar period value was obtained with high confidence 
from later Swift and XMM-Newton observations, and it also 
coincides with the rise in brightness and the appearance of an ap¬ 
parent emission line spectrum with many unknown line features. 

Since after day ~ 29, the 35s period was always associated 
with SSS emission, this coincidence may indicate that the 
strange soft component with emission lines (Fig. [T]i is also 
associated with SSS emission, even though the spectrum does 
not appear to be a typical SSS continuum spectrum. In this 
context, it is noteworthy that the low-resolution simultaneou s 
EPIC spectrum can be reproduced by a blackbody (iNes ilMl. 


A Chandra observation taken on day 39.7 after outburst 
(top panel of Fig. HJi contains a peak at 37.5 s that proved to be 
insignificant. This observation was taken during an episode of 
low c ount rate during an early variability phase (lOsborne et al.l 

IMih . 


During the XMM-Newton observation taken on day 54, the 
peak of the overall X-ray brightness evolution dOsborne et al.l 
1201 ih . the 35s period was detected at a significance > 99.9%. 
The extreme brightness of RS Oph at this time led to difficulties 
in telemetering the data to the ground, and the EPIC data can 
only be used for spectroscopy. The RGSl light curve could 
not be used for the same reason; only the RGS2, which had a 
lower count rate owing to a failure of a CCD chip (that would 
have recorded the dispersed photons from the peak of the SSS 
spectrum), contained a controllable number of telemetry drop 
outs. In Fig. |2] we show the time evolution of the period (left) 
and of the amplitude (right) from the RGS2 data. The signal 
power varies considerably, as shown by the two power spectra 
in the top panel, extracted from the time intervals marked in the 
right panel of the left plot. The orange open histogram contains 
two much weaker signals at ~ 34 s and ~ 36 s. In the right panel, 
the evolution of the phase-folded light curves is shown on the 
same (vertical) time axis as the left plot. Since the start and stop 
of the time segments do not coincide with the same epochs, 
we have artificially moved each curve such that the peak is 
always at phase 0.25. The right plot in the right panel shows the 
evolution of the amplitude. 
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period (s) 



phase 


Fig. 2. RS Oph on day 54. Left: Time map of the 35s period in the XMM-Newton RGS2 light curve, see description of right panel 
in Fig. [1] Right: Time map of folded light curves, where each curve was artificially moved in phase for the peak to occur at phase 
0.25. In the top panel, an example is shown with the sine curve fit from which the modulation amplitudes with uncertainties were 
determined. In the right panel, the evolution of amplitudes is shown. 


We studied the relations between the period, power, ampli¬ 
tude, and the mean count rate in 38 time segments of the full 
light-curve, excluding six segments in which the amplitude only 
yields an upper limit. In Fig. [2 the values of these parameters 
are plotted against each other. In the left column, the three 
period parameters are plotted against the mean count rate. In 
the top row, the 95% (2cr) detection threshold (corrected for 
oversampling) is indicated by the horizontal line. The signal 
power does not seem to depend on the count rate. Out of the 
32 time segments with positive amplitudes, 30 yield > 2cr 
detections of the 35s period, and these are used in the two panels 
below. The 95% error bars included for the periods in the middle 
left p anel are twice t he standard deviation, calculated with Eq. 
(4) in|L arssoni (Il996h . The error-weighted mean period is given 
with the horizontal line with the value and standard deviation 
given above this line. The errors are smaller than the widths of 
the peaks, which are driven by windowing and thus the duration 
of the time segments T. The frequency resolution scales with 
l/T^, and at the period of 35s, the expected width of the peak is 
dp = IT ~ Is. The variations of the period are slightly larger 
than the 95% error bars and thus may be real, although we note 
that the large uncertainties in amplitudes have not been ac¬ 
counted for when calculating the period errors. Furthermore, the 
distribution of periods, shown in the small inset, indicates that 
the variations are normally distributed, because the distribution 
can be fit by a Gaussian. The distribution and formal errors do 
not provide information to conclud e whether or not the pe riod 
varies with time, as concluded by lOdendaal et al.l (1201 4l) for 
the 67 s period in Cal 83. In the top panel of Fig. |2 the power 
spectrum derived from the complete light curve is shown by the 
thin black line, and the higher period resolution achieved by the 
longer baseline of the total light curve reveals multiple peaks. 
This would support that the period is not constant in time. 


In the bottom left panel of Fig. [2 the amplitudes in units 
of RGS2 count rate are shown against the average count rates 
for each time segment. The dashed and dotted lines indicate 
the lowest and highest amplitudes relative to the respective 
count rate in per cent. The amplitudes range between 2.3-11 
counts per second (1.7-7% of the total count rate) without any 
systematic trend with count rate. Uncertainties in amplitude 


range between 30% and 90%. 


In the right column of Fig. [2 relations between the am¬ 
plitude, power, and period can be studied. Obviously, small 
amplitudes are more difficult to detect because they yield low 
detection powers. Some outliers with low signal power but 
high amplitude can be caused by a high degree of additional 
non-periodic variability, which in a folded light curve could 
result in a large amplitude. The uncertainty in amplitude is 
marked with the horizonal grey line in the top left panel and is 
larger than the deviations from the main trend. In the middle 
right panel, the variations of period with amplitude are shown 
and it seems that the larger variations in period coincide with 
low amplitudes. This trend seems even stronger in the bottom 
right p anel, which plots power versus period. lOsborne et al.l 
(1201 Ih showed that the period stability increased with time 
through the SSS phase. 


Another Chandra observation taken on day 67 after outburst 
may present some evidence for the 35s period at the 1-cr level 
(see Fig. HJi. The colour scheme and thresholds for testing 
include the correction for oversampling and for multiple testing, 
although the latter is conservative and does not recognise 
that this p articular perio d was clearly detected in other ob¬ 
servations. lOsborne et alJ (1201 ih reported that the 35s period 
disappeared on day 58.8, about a week before this observation 
was taken. If the period is real, then it might be related to a short 
re-brightening tha t seems evident from the Swift light curve 
shown in Fig. 2 in lOsborne et alJ (1201 ll) . On the other hand, the 
Chandra observation contains no re-brightening, in particular 
not around the time of this marginal detection. 


Only one 4400-sec exposure was taken with the XMM- 
Newton optical monitor on day 54 with the UVW2 filter. While 
the exposure was taken in fast mode, the source was not cen¬ 
tred within the small readout window, leaving most of the PSF 
outside of it. The light extracted from the remaining part of the 
source contains no significant periodic signal. 
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Fig. 3. RS Oph on day 54: Relations between power (top row) and period (middle row) versus mean RGS2 count rate (left column) 
and modulation amplitude (right column). The amplitude is mean to peak in units of RGS2 counts per second. A typical uncertainty 
range in amplitude, derived while fitting sine curves to the phased light curves, is shown with the grey horizontal line in the left 
corner. The 95% (2cr) detection level (corrected for oversampling) is marked by the horizontal line in the top panels. In the middle 
panels, the error-weighted median period with standard deviation of the peri ods is marked w ith the horizontal solid and dotted lines, 
respectively. The 95% error bars on periods were calculated with Eq. (4) in lLarssonI (Il996h (see text). We show 30 out of 38 data 
points for which the signal power exceeds a 95% detection and a sine fit yields a positive amplitude. The distribution of periods is 
shown in the small inset where a Gaussian is added for comparison. In the middle right panel, period versus amplitudes are shown, 
suggesting that larger variations in period occur with smaller amplitudes. In the bottom left panel, amplitudes are shown versus 
count rate, demonstrating that there is no correlation. Ranges of amplitudes relative to count rate are marked by the dashed and 
dotted lines. In the bottom right panel, periods versus power are shown, which also seem to show that larger variations in period 
occur at lower values of power. 


4.2. KTEri 

KTEri was the second nova in which a transient 35s period was 
seen in Swift monitori ng observations during the SSS phase 
dBeardmore et al.ir2010t) . Eour Chandra observations were taken 
during the eruption on days 71.3, 79.3, 84.6, and 158.8, and we 
searched the light curves for periodic oscillations. Owing to the 
high brightness of the source and with the main goal being to 
obtain high-resolution grating spectra, the observing times were 
rather short. In such short exposures, transient periodicities can 
be missed easily. 

We detected the 35s period only during the early and late 
parts of the observation on day 159 (see Eig. |5]l. The width of 
the peak is consistent with the highest possible resolution power 
in period, associated with the selected duration of each time seg¬ 
ment. Extracting the power spectrum from the total light curve 


yields a higher resolution in period, as can be seen from the thin 
black line in the top panel. The double nature of the peak sug¬ 
gests that the period was variable. 

4.3. V339Del 

The third time a transient short-period signal was seen in Swift 
observati ons of a nova wa s V339 Del (on days 77.5 to 88.6), 
in which iBeardmore et al.l (1201 3h reported a 54s period; this 
was conhrmed to have also be en present i n an X MM-Newton 
observation taken on day 99 (iNess et al.l 12013bh . While the 
XMM-Newton EPIC-pn camera suffered a full scientific buffer 
at almost all times and could be used only for spectroscopy, 
simultaneous light curves from the MOS and RGS instruments 
provide independent evidence. The time evolution of power 
spectra derived from the MOS2 data is illustrated in Eig. 
where we can see that the signal was not constant in power. 
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period (s) 

Fig. 4. RS Oph on day 67; Time map around the 35s period in 
the Chandra HRC light curve, see description in the right panel 
of Fig. [U The strongest peak coincides with the 35s period seen 
in earlier Swift and XMM-Newton observations (see Fig. 13. The 
power (conservatively corrected for multiple testing and over- 
sampling, see Appendix) only yields a 1-cr detection. 



Fig. 5. KT Eri on day 159: Time map around the 35s period in the 
Chandra HRC light curve, see description in the right panel of 
Fig.m The thin solid line in the top panel is the power spectrum 
from the entire observation and shows two > 99.9% significant 
peaks. This indicates that the period is not constant in time. 


similar to RS Oph. We also analysed the two Chandra obser¬ 
vations taken on days 87.2 and 114 and found similar results. 
We performed the same parameter studies as for RS Oph for 
all time segments from all three observations and show the 
results in Fig. |7] The XMM-Newton and Chandra data can 
clearly be distinguished by the lower Chandra count rates - 
owing to the smaller effective areas - and consequently smaller 
amplitudes, which are given in units of count rate. The mini¬ 
mum and maximum ratios of amplitudes versus count rates are 
indicated in the bottom left panel with dashed and dotted lines, 
respectively. The relative amplitudes are consistent for Chandra 
and XMM-Newton data. In V339Del, as in RSOph, no clear 
relation between the power and the brightness of the source can 
be found. The variations in period are slightly larger than the 
95% uncertainties in period and thus may likewise be real. The 
modulation amplitude of signals detected with more than 95% 
confidence varies between 3.1-10% of the count rate, clearly 
yielding more significant detections (higher powers) with higher 
amplitudes (top right panel of Fig. |7]). The slope in the top right 


panel is steeper for the Chandra than the XMM-Newton data 
owing to the difference in sensitivity and background. Again, 
like RS Oph, the amplitude variations are not related to the 
brightness of the source (bottom left panel). The bottom right 
panel suggests that the largest scatter in period is found for the 
lowest power. 

The XMM-Newton optical monitor was operated in fast 
mode, taking ten exposures with the UVW2 filter. The count rate 
ranged around 140+lOs ', consistent with Poissonian noise, and 
we found no periodic oscillations. 

4 . 4 . LMC 2009a 

We found a new system with a 33s period in two out of four 
XMM -Newton observations of nova LMC 2009a (iNess et al.l 
l2014t) . In Fig. |8] the time maps are shown for all four observa¬ 
tions. In the first observation, taken on day 90 after outburst (top 
left panel), no periodic signal is seen. In the second observation, 
taken on day 165 after outburst (top right panel), a significant 
signal with a period of 33s can be seen at ~ 8 hours of elapsed 
time. The same period is more clearly detected in the third 
observation taken on day 197 (bottom left panel). In the last 
observation, none of the time segments contains a significant 
signal. We note the single peak in the power spectra derived 
from the respective full light-curves. While a peak at 33s is ap¬ 
parent, it is not statistically significant. In the bottom left panel, 
the total power spectrum contains multiple peaks around 33s, 
clearly indicating strong variability of the period. We performed 
the same studies on correlations as those shown in Figs. |3] 
and |7] using 22 data points with > 2cr detections and positive 
amplitudes and find similar results of significant variability 
in the period, while variations in power and amplitude do not 
correlate with the count rate. 

As can be seen from Fig. [8] the fraction of time that the 
33s period is detectable in LMC 2009a is much lower than for 
the periods in RSOph, KTEri, and V339Del (see Sect. 14.7t . 
Perhaps this explains why it is not found in the SwiftfXWV data 
(Beardmore in prep); Swift might have observed only during the 
times when the period was absent. 

Each observation was accompanied by 20-30 fast mode ex¬ 
posures with the optical monitor, all using the UVW1 filter. The 
count rates decreased with time starting from 6.3 +1.3 s“' on day 
90.4 through 2.4 + 0.7s“' on days 165 and 196 to 2.1 ±0.6s“' on 
day 230. Lrom each light curve we constructed the power spec¬ 
trum, but found no significant signal, even when combining all 
exposures into a single light-curve. 

4 . 5 . Cal 83 

Lor the persistent SSS Cal 83, lOdendaal et akl (12014l) have re¬ 
cently reported a 67s periodic modulation. We have produced 
our pe riod time maps for 19 XMM-Newton observations, for 
which lOdendaal et al.l (l2014ll report that the source was not 
in the off state, and present an example in Lig. 0 for ObsID 
05 06531701 that yields the same result as that shown in Lig. 8 
of lOdendaal et al.r(l2014l) . 

We repeated the same correlation studies of the power of 
the period and period with count rate that we have presented 
for RSOph and V339Del. Lor Cal83, we have a much larger 
sample, which is shown in Lig. (TO] In the top panel, we include 
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Fig. 6. V339 Del on day 99: Same as Fig. |2] focusing on the 54s period in V339 Del based on the XMM-Newton MOS2 light curve. 
The simultaneous combined RGS1 and RGS2 light curve yields a consistent result. 



count rate (s ') power 

Fig. 7. Same as Fig. [3] for the Chandra and XMM-NewtonlMOSl observations of V339Del. We show 94 out 266 data points in 
the middle and bottom rows that yield > 2cr detections and positive amplitudes. As in RS Oph (Fig.O, no relation between period 
properties and count rate is apparent, and the period seems to be variable. The modulation amplitude clearly correlates with the 
power of the period because modulations with a larger amplitude are easier to detect. The Chandra and XMM-Newton data yield 
different slopes for amplitudes in unit counts per second, while for count rates, a single slope is encountered. In the bottom left panel, 
we show that the data from both telescopes yield similar ranges of relative amplitudes. The bottom right panels appear similar to 
RS Oph (Fig. [3]l, with largest variations in period occurring with low values of power. 


the results from all time segments. The variations in period 
are again slightly larger than the 95% errors and thus might 
be real. At intermediate count rates, the period seems to be 


anti-correlated, but this would be based on too few data points 
(middle row). We also searched for periodic variations of the 
power on amplitude on longer time scales, but found no clear 
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Fig. 8. LMC 2009a on days 90, 165, 197, and 230: Time maps around the 33s period in four XMM-NewtonIpn light curves, see 
description in the right panel of Fig. [1] The vertical (time) axes are different in each plot; they only cover the respective exposure 
times. 



period (s) 

Fig. 9. Cal 83 on 2009-05-30.4: Time map around the 67s period 
in one of 19 XMM-NewtonIpn light curves, see description in the 
right panel of Fig. [T] The time evoluti on is consistent with th at 
shown in a different format in Fig. 8 in lOdendaal et al.l (l2014l) . 


patterns in the power spectra of signal powers and amplitudes. 

We checked all available XMM-Newton observations with 
the optical monitor, which contain exposures in U, B, UVWl, 
UVM2, and UVW2 hlters, but found no periodic signal in any 
of these exposures nor in combined light curves for each hlter. 


4 . 6 . Search for other candidates 

We screened all power spectra from the datasets listed in Table[T] 
for any signihcant signals that might account for multiple testing 
and oversampling. With these corrections any peak at any period 
that rises above a given threshold is a strong candidate. We 
searched for signals that rise above the 95% threshold, but found 
only the hve sources already discussed above, thus no other 
candidate can be identihed. 


4 . 7 . Fraction of observing time for which short-period 
osciiiations are detectabie 

To compute the fraction of time for which the signals are 
detectable above a given threshold, we recomputed power 
spectra for adjacent, non-overlapping time segments of 1000s 
duration. We hrst dehned a detection threshold and summed 
the exposure times of all time segments in which a period is 
detected above the given threshold. We then computed the ratio 
of the summed exposure time by the total available exposure 
time. The resulting fractions depend on the chosen signihcance 
level, yielding higher fractions for lower detection thresholds. 

For each target, we probed a range of detection levels, and 
in Fig.[TT]we show fractions of times as a function of detection 
threshold for the hve sources with detected periodic signals. 
Since the period was not necessarily present in all observations, 
the global ratio strongly depends on the coverage of the total 
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Fig. 10. Same as Fig. [3] with data collected from 20 observations of Cal 83, divided into a total of 578 1000s time segments. 17 
of these intervals yield zero count rate and were excluded. Furthermore, all data points yielding only upper limits in amplitude 
were excluded, which left 202 data points, of which 34 time intervals yield a > 95% detection of a periodic signal. The expected 
correlation between modulation amplitude and signal power (top right) is contaminated by a large number of overestimated ampli¬ 
tudes that are most likely caused by additional, non-periodic variability. These deviations are consistent with the large uncertainties 
in amplitude. The 95% threshold is corrected for oversampling, but not for multiple testing. 31 of 34 data points that are found 
above this threshold a nd yield a successful sine fit are used in the lower two rows. The variations in period of about +3s reported by 
lOdendaal et al.l (l2014l) can be recognised in the middle row and bottom right, but do not seem to correlate with count rate, amplitude, 
or signal power. 


evolution, which is generally much better determined with 
Swift. To study the short-term behaviour, we also include in 
Fig. HU the ratio for the observation with the highest fraction 
with dotted lines and labels in italics for each source. For all 
sources, solid lines indicate that the transient signal is episodic. 
These estimates can be refined with the Swift data, which will 
be presented by Beardmore et al. (in preparation). 


4.8. Search for patterns in the transient signal 

The short-period oscillations are clearly transient in all five sys¬ 
tems with detected periods, and we investigated whether there 
are any patterns in the power of the signal. For example, pe¬ 
riodic variations could be caused by the rotation of the white 
dwarf if the oscillations originated in an accretion spot. For each 
source, we created a time series of the signal power and period 
value from all time segments. We computed a power spectrum 
based on the central times of each time segment and the associ¬ 


ated parameters. None of the resulting power spectra reveal any 
significant signatures. 

5. Discussion 

We have four confirmed novae with periodic or quasi-periodic 
transient signals ranging between 33 and 54 s and one persis¬ 
tent SSS with a periodic signal of 67 s. For these five systems, 
we list in Table [2 the periods, limiting count rates, and rela¬ 
tive amplitudes of the modulation for 2-cr detections, the opti¬ 
cal speed class parameters f 2 and t^, and fractions of time when 
detected. The limiting count rates and amplitudes are the low¬ 
est value from the sample of time intervals when a period was 
detected with more than 2cr significance. Although this sample 
is too small for searches of systematic trends, it covers a wider 
range of parameters than the time series of individual systems il¬ 
lustrated in Figs. [2111 andfTOl No systematic trends can be seen 
between the period and the count rate or amplitude, but the nova 
speed class seems to scale with the period, yielding longer peri¬ 
ods for slower novae. The speed class is known to scale with the 
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Fig. 11. Fractions of observing time for which the strongest peri¬ 
ods are detectable as a function of detection threshold. The solid 
lines are the average fractions of time derived from all available 
observations, while the dotted lines (with labels in italics) are the 
fractions of time derived from the single observation in which 
the period is seen most frequently. For this plot, power spectra 
from directly adjacent, non-overlapping time intervals have been 
used. For no system, the period was detected 100% of the time, 
showing the transient nature. Furthermore, the signal is episodic 
because the fractions are not the same in all observations. 


mass of the underly ing white dw arf, yielding higher masses for 
shorter speed class (lLiviolll99^ . and shorter periods thus may 
be present in systems with more massive white dwarfs. 

5.1. Detectability of the period during Thomson scattering? 

In this section we show that Thomson scattering does not reduce 
the detectability of the period. Recent studies of the eclipsing 
recurrent nova U Sco have revealed that about 50% of the con¬ 
tinuum emission from the white dwarf is still visible during total 
eclipse. The shape of the continuum spectrum seems to have 
been preserved, while emission lines dominate the spectrum. 
This indicates substantial effects of Thomson scatte ring of the 
contin uum while the emission lines are not eclipsed dNess et al.l 
I2OI2I1 . A few other SSS spectra show similar characteristics 
of reduced but d etectable continuum emission with emission 
lines on top, and iNess et aTI (l2013ah have classified this type 
of SSS spectrum as the SSe subclass: ’e’ for emission lines as 
opposed to SSa with continuum-dominated spectra in which 
deep absorption lines can be seen. 

IChang & Kvlafisl (Il983h have studied the effects of scatter¬ 
ing of photons on the spectra and variability patterns of X- 
ray sources. Since Thomson scattering is an energy-independent 
process, their calculations of spectral changes are not of inter¬ 
est to us, and we focus on their predictions of reductions of the 
amplitude of oscillations caused by scattering. These reductions 
are shown in their Figs. 7 and 8 and depend on the ratio of the 
period p to the light travel-time, p c/ro and the optical depth, 
T = Jq HeCTfdr with cTf = 0.6652457x 10“^"^ cm^ the cross-section 


for Thomson scattering by electrons, and rig the electron density. 
For our purposes, we focus on the models for a uniform sphere 
because we most likely examine a continuous outflow and not 
an abrupt ejection event that would produce a thin shell. We 
have observed amplitudes of about 10% of the count rate and 
can thus concentrate on the rang e of amplitudes v i ewed from 
infinity log(Aoo) > -1 in Fig. 8 in IChang & Kvlafisl (Il983h and 
thus log(p/f) > 0.5 with f = ro/c the light travel-time. While 
this applies only to the case t = 10, lower values of r yield 
some 0.2 dex lower value s of log( p/t) for log(Aoo) > -1, which 
can be seen from Fig. 7 in IChang & Kvlafisl (I l983h . We thus fo¬ 
cus on the range log(/7/f)^0 and thus < c p. For the period 
range 33-67 s, this constrains the size of the scattering region to 
ro < (14-29)R0. 

Depending on the assumed optical depth t, this allows us to 
place a constraint on the electron density tig of this plasma. With 
T = tigCTtr for a uniform sphere and r < c p, the constraint on 
density is thus rig > T/(cr, c p). 

In Fig. [12] the derived densities for a range of optical depths 
T =1-5 (assuming that 1-37% of photons escape unscattered) 
are illustrated for a period range of 25-85 s. In the extreme 
case of an intrinsic amplitude of 100%, Thomson scattering 
can reduce the amplitude to 10% of the count rate only if the 
electron density of the scattering medium exceeds 10'^ cm“^. If 
the intrinsic amplitude is much below 100%, then the required 
densities would have to be even higher. A known intrinsic 
pulsed fraction and optical depth t would allow us to work out 
both the size and density of the scattering region by the signal 
contrast degradation, but the intrinsic amplitude is probably not 
much higher than the observed one. 

We thus conclude that Thomson scattering does not effec¬ 
tively reduce the detectability of the periodic signal, and in fact 
we also found a periodic signal in the SSe LMC 2009a, even 
though it is the system with the smallest fraction of total ob¬ 
served time with significant detections and the smallest ampli¬ 
tude. 


5.2. Rotation period of the white dwarf 

The long-term pres e nce o f the 67s X-ray period in Cal 83 has 
led lOdendaal et al.l (l2014l) to consider the possibility that the 
rotation period of the white dwarf is the main driver of the 
oscillations. While typical rotation periods of white dwarfs 
are much longer, the fastest rotation period in a white dwarf 
binary, RXJ0648.0-4 418/HD49798, was observed to be 13.2s 
(IWarner e t al.||2003b|), and a rotation period of 67s is thus pos- 
sible. lOdendaal et^ (l2Q14ll argued that a long mass accretion 
history allows the white dwarf to be sufficiently spun up. The 
greater challenge would be to explain the observed changes 
in period by +3s. Because of the large amount of inertia of a 
white dwarf, this drift cannot be the result of accretion-induced 
spin-up or -down. In addition, there i s no relation with the 
orbital period o f 1.0475 + 0.00004 days (ISchmidtke et alJl2004l : 
lOdendaal et al.llTOTih . which excludes Doppler shifts imposed 
by the orbital motion of the X-ray emitting plasma, for example. 
Furthermore, the possibility was discussed that the observed 
period is the beat between the rotation period of the white dwarf 
and the Keplerian period of the denser X-ray emitting plasma, 
but the required rotation period of the white dwarf would then 
have to be even shorter, with an estimated value of 4-12 s, which 
would set a new record, althou gh this is still longer t han the 
break-up period of 3.4 s. Finally, Odendaal et al.l (1201 4h argued 
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Table 3. Properties of systems with known periodic signals 


Target 

period® 

erf 

lim 

A“ 

A“ 

h 

h 

duty cycle”' 


(s) 

(s-‘) 

(s-') 

(%) 

(days) 

(days) 

(%) 

EMC 2009a 

33.2 ±0.1 

45 

1.5-2.9 

2.9^.8 

4(1) 

11(4) 

11.3/38.1 

RS Oph 

34.9 ±0.2 

148 

2.3-11 

1.5-7.0 

7.9®/6.2<5) 

17.1® 

38.8/74.5 

KTEri 

35.0 ±0.1 

70 

1.5-4.2 

2.5-3.9 

6.6®/6.2<'^’ 

14.3® 

16.9/61.4 

V339Del 

54.1 ±0.3 

176 

6.6-13.8 

3.4-7.6 

12 ± 2® 

23.5® 

49.2/67.8 

Cal 83 

66.8 ± 0.5 

4.0 

0.6-1.0 

7.0-15.7 

- 

- 

18.7/61.8 


“Errors on periods are the standard deviation. 

^Minimum count rate during time intervals when periods with > 2cr are found 

“Amplitude range in counts per second (Col . 4) and in % of count rate (Col. 5) that yields > 2cr detections _ 

t 2 and ts: from (lllSc hwarz et"^ ll201 ih . (2l lHounsell et al.l ll2010h . (.Sl lMunari et all ll2013h . (4l lMason & Walt^ ll2014h . (.Sl lMunari et al.l ll2007h . 
(61 llmamura & Tanab?1^12ll 

“^Fraction of observation time (all observations combined/in the single observation with highest fraction) detected at > 95% (2<r) 



Period [s] 


Fig. 12. In the presence of Thomson scattering, the observable 
amplitude of a signal will be lower than the original modula¬ 
tion amplitude, depending on optical depth, t, and density, n. 
We show estimated densities below which periodic oscillations 
would lead to reductions of the amplitude (and thus detectabil¬ 
ity) by 64-99%, depending on three values of optical depth 
T =5-1 that would lead to 1-37% of photons escaping unscat¬ 
tered; see Sect. O 

that the X-ray emitting plasma must originate from an extended 
envelope, some two to three times the radius of the white dwarf, 
and if that envelope does not rotate synchronously, the spread in 
the 67 s oscillations can result. However, this is not quantified, 
and the cause of the modulation within this gas remains unclear. 

Within the sample of detected periods, the 67s period in 
Cal 83 is the longest one, and interpretating this as the rotation 
period would imply that the white dwarfs in the Galactic novae 
in our sample would be even shorter. Spinning up to such short 
rotation rates impli es even longer mass-accretion histories. 
lOsborne et al.l (1201 ih also argued that the period changes are not 
consistent with realistic torques on the white dwarf. The changes 
in period by +3 s argue against rotation of the white dwarf, al¬ 
though the argument about non-synchro nous modulation within 
a hydrogen gas envelope presented by lOdendaal et al.l (l2014t) 
cannot be fully discarded. Thi s non-synchron ous region may 
be similar to that proposed by IWarnerl (l2004h . who explained 
’dwarf nova oscillations’ by appealing to weak m agnetic fields 
that incompletely enforce co-rotation; see also IWoudt et al.l 


(l2010l) . However, it remains unclear how the rotation period 
would manifest itself as an observable period, for instance, 
whether it might be caused by periodic absorption patterns. 

iBeardmore et al.l (120101) noted that the presence of the same 
periods in KT Eri and RS Oph would appear to argue against a 
rotation-based origin for this modulation. Meanwhile, a possible 
residual nuclear-burning white dwarf pulsation might be more 
constrained in frequency, which would also explain the small 
range of detected periods in our sample. 

The rotation of the white dwarf might be an interpretation 
for the observed oscillations if some speculative arguments were 
accepted, but we cannot identify under which circumstances this 
period can be seen in X-ray observations. In particular we can¬ 
not explain why the signal is transient and why it is seen in 
such a small fraction of systems observed. If magnetic fields are 
needed, as in intermediate polars, the question arises why these 
oscillations are seen only during the SSS phase, while interme¬ 
diate polars show the rotation period in optic al and X-ray light 
curves in quiescence, for example, V4743 Sgr (iDobrotka & NessI 
l2010h . An independent measurement of the rotation period for 
one of our five systems with short-period oscillations would ro¬ 
bustly confirm or reject the possibility that the rotation period is 
the main driver for these oscillations. Meanwhile, other possibil¬ 
ities such as pulsations exist. 

5.3. Pulsation modes 

An attractive explanation for the observed oscillations is that it 
is a surface ^-mode that is triggered by the e-mechanism. This 
possibility gains some favour from considering the dynamical 
time-scale of the white dwarf atmosphere. In the e-mechanism, 
the contraction and compression of a mode leads to changes 
in the nuclear burning rate. If this energy input exceeds the 
cooling during the subsequent rarefaction of the mode, it will 
be unstably excited. Since the oscillations are preferentially 
observed in the SSS following nova outbursts, when nuclear 
burning may be occurring at the surface, the conditions may be 
met for this to occur. 

The e-mechanism has been considered for g-modes in the 
context of white dwarfs following the production of planetary 
nebulae (lKawalerlll988l) . as well as for neutron stars that are 
stably burni ng their accreted fuel si milarly to the way this occurs 
in an SSS dPiro & BildstenI l2004h . Since this scenario is for 
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physical conditions somewhat different from the nova systems 
studied here, it is helpful to provide some basic arguments about 
what types of periods we expect for surface g-modes. This 
also demonstrates that the mode periods roughly depend on 
the properties of the white dwarfs and thus indicates that these 
modes may be useful probes of their structure. 


When the white dwarf is in an SSS state, there is an abrupt 
change in the internal profile (in composition, temperature, and 
density) where hydrogen is being burned. There is an associated 
strong buoyant force at this depth, so that we can consider modes 
that are riding along, unable to penetrate this region. If the burn¬ 
ing occurs at a depth H, which is much much lower than the 
radius of the white dwarf R, this mode’s di spersion rel ation will 
be similar to that of a shallow water wave (lFaberlll9^^ 

— , ( 1 ) 

P 

where at is the mode frequency, g is the local gravitational 
acceleration, g » GMjR^ when H R, = /(/ -i- l)/R^ is 
the transverse wave number squared, where I is the spherical 
harmonics quantum number, and Ap/p is the fractional change 
in density at the burning depth. For simplicity, we assume 
Ap/p ~ 1 in our discussion below, which corresponds to a 
maximally discontinuous jump. 


While g and k depend on the mass and radius of the white 
dwarf, H depends on the thermal profile of the surface layers. 
The thermal profile of the white dwarf is set via flux-limited dif¬ 
fusion, which obeys the equation 

L ^ 4acT^ dT 
AnR^ 3 kp dr ’ 


where L is the luminosity, a is the radiation constant, k is the 
opacity, T is the temperature, and p is the density. We assume k - 
0.34 cm^g^' is a constant as set by electron scattering for solar 
composition material, which is roughly correct for the hot SSS 
state. Using hydrostatic balance, dPIdr - —pg, this expression 
can be rewritten as 

L ^ AacgT^ dT 

AnR^ 3k dP' ^ 

Assuming that L is constant with depth (which is a good ap¬ 
proximation above the burning depth in the SSS), this can be 
integrated to yield 


r(P) = 



(4) 


where Ledd = AttGMcIk is the Eddington luminosity. The layer 
thickness is then given by 

) (—) , (5) 

pmpg pmpg\aj \LEdd/ 

where k^ is the Boltzmann constant, p is the mean molecular 
weight, nip is the proton mass, and we assume an ideal gas 
equation of state. 


via theoretical calculations. For now, we know that the pressure 
must be lower than the pres sure at ignition of the layer (for exam¬ 
ple, see the calculations of iTruran & Liviolll986h . so for demon- 
strational purposes we set P « 10'^ dyn cm“^, so that 


H ~2xl0' 


8 


v2 X lO^cms 


-2 


r( 


lO'^ync 


1/4 



cm, (6) 


where we use g as appropriate for a roughly a; 1 .3Mq, similar to 
the kinds of masses inferred for these systems. Substituting this 
into the shallow water wave dispersion relation, the estimated 
period is 


/ p \ 1/8 

mode y 1018 dyn cm“2 ) 




(7) 


where we take / = 1 since this is the lowest order (and thus eas¬ 
iest to observe) angular mode, and Ap/p ~ 1 due to the strong 
compositional stratification following the burst. This result high¬ 
lights that the strongest dependency of Pmode is on R and that 
for typical parameters associated with recurrent novae the pe¬ 
riod is too low. The two possible solutions that could reconcile 
this period with the observed periods are (1) the white dwarfs 
are much less massive and have a larger radius, or (2) the higher 
order radial modes are excited instead of the lowest order mode 
we consider here. The first possibility is unlikely given the con¬ 
straints on the white dwarf masse s from the recurrence times and 
nova durations (IWolf et al.ll2013h . The second possibility would 
give a better match in period since the higher order radial modes 
increase with the number of radial nodes. This would then re¬ 
quire some explanation for why the higher order mode is excited, 
which requires a more detailed theoretical investigation. 

If this mode explanation for the oscillations is confirmed, 
then the period of these modes provides a constraint on the 
mass of the white dwarf. This could be useful for determining 
whether these nova systems are more massive than typical white 
dwarfs, which may indicate that they are growing from accretion 
over time. Signs of such mass growth may have implications for 
whether some of these systems are attractive as Type la super¬ 
nova progenitors. 


6. Summary and conclusions 

We focused on analysing continuous light curves obtained with 
XMM-Newton and Chandra. We constructed power spectra that 
probe the period interval 25-lOOs and tested 500 frequencies 
on an equidistant frequency grid, which corresponds to an 
oversampling of 16. We estimated detection probabilities (ap¬ 
pendix) based on simulations of 10,000 synthetic light curves 
with random Poisson noise. We furthermore broke each light 
curve into 1000s time segments to study time variability of any 
periodic signal. We corrected the detection probabilities for 
multiple testing by the conservative Bonferroni correction. 

We confirm the transient nature of the 35s periods in RS Oph 
and KTEri and the 54s period in V339Del and determined the 
fractions of time that the respective periods are detected at 2cr 
significance as listed in Tabled The 67s period in Cal83 and 
the 33s period in EMC 2009a are detected less frequently. 


The thickness PI therefore depends on the luminosity L, 
which can be inferred from observations, and the pressure P at 
the burning depth. This pressure should be related to L and the 
requirement of steady burning. This therefore can be addressed 


Our dynamic power spectra illustrate variations in power 
and amplitude on time scales of fractions of an hour. We also 
found multiple peaks in the global power spectra for each 
observation, suggesting variations of the period, although the 
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uncertainties in periods are too large to confirm that. None of 
these variations correlate with the X-ray brightness. 

The amplitudes are found to vary by about <15%. We 
checked the possibility that the amplitudes (and thus detectabil¬ 
ity) may be degraded by Thomson scattering, which could 
explain why so few sources are found with a period and why 
it is transient. We found that the density of the scattering 
medium would have to exceed lO'^cm^^. It depends on the 
unknown optical depth and the intrinsic amplitude. The intrinsic 
amplitude is probably much lower than 100%, thus requiring 
even higher densities, so we conclude that Thomson scattering 
has no significant effect. 

Potential origins of the observed oscillations are the rota¬ 
tion of the white dwarf or intrinsic pulsations. Only speculative 
arguments might explain the observations with rotation of the 
white dwarf, and we thus focused on pulsations. While periods 
in the 10-100 second range are predicted by nonradial oscilla¬ 
tor modes or pulsations in isolated white dwarfs, the validity of 
these models under conditions of active nuclear burning during a 
nova outburst is difficult to estimate, and new models need to be 
calculated to account for the higher luminosity and the physical 
conditions of nuclear burning near the surface. 
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Appendix A: Appendix 

The period searches were performed using power sp ectra calcu¬ 
lated with the method of [Horne & BaliuiiM (Il986h . where the 
power is norm alised to the total variance of the entire data set. 
IScarglel(ll982h emphasised that surprisingly large spurious peaks 
can occur in a power spectrum and investigated the probability 
that a signal with a given power z is a spurious signal (false- 
alarm probability). His Eq. (14) computes the probability for 
a given power peak to be due to a chance noise fiuctuation, 
including the statistical penalty for inspecting a large number 
of frequencies, Nf. Based on this definition of the false-alarm 
probability, P^p, we define the reverse, a detection likelihood 
-C = 1 - Pfap, thus; 

X = (l-e-^)^<, (A.l) 

from which the power z that yields a given detection likelihood 
-C can then be inferred from 

z=-/n(l(A.2) 

For example, in a power spectrum consisting of /Vf = 500 tested 
frequencies, the required normalised power that yields a 99.9% 
detection probability is z = 13.2. However, in the discussion that 
follows we determine the most appropriate choice for Nj, and 
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Fig.A.l. Simulations based on N - 10,000 synthetic light 
curves of 1000s duration with Poisson noise around various 
given mean count rates to estimate the likelihood of random de¬ 
tections at four levels of confidence. Power spectra of the ran¬ 
dom light curves were calculated over the period range 25-100 s, 
testing 500 frequencies. Left panel: Example of 20 counts per 
second as mean count rate. We show the maximum powers in the 
order of trial number in grey and sorted by power value with the 
thick black line. The sorted curve was used to extract the power 
thresholds given in the legend. In the right panel we show the 
results for the four confidence levels as a function of mean count 
rate with connected black bullet points in comparison to the the¬ 
oretical thresholds according to Eq. (IA.21 l. with Nf = 500, in¬ 
cluded with the dotted blue lines. 


hence the detection threshold, in an oversampled periodogram. 

In each 1000s time segment studied in this article, the 
period range 25-100 s is probed, corresponding to the frequency 
interval 0.01-0.04Hz and thus a frequency range df - 0.03 Hz. 
The sampling above from which no further information can 
be gained is to test at the Eourier frequencies, thus with a 
frequency bin size of 1/f, d/ x f = 30, within the 0.03Hz 
frequency range. Meanwhile, with Nf = 500 frequencies, we 
oversample with a factor 16. By oversampling, we are more 
sensitive to periodicities that might fall midway between the 
Eourier frequencies. 

As a test of the detection thresholds defined above, we 
performed simulations by generating synthetic light curves with 
random Poisson noise around a range of given mean count 
rates. Eor each given mean count rate, we generated a sample 
of 10,000 synthetic light curves of 1000s duration, the same 
duration as was chosen for the time segments into which we split 
the observations. We then computed power spectra under the 
same conditions as for the data, thus probing a period range of 
25-lOOs with Nf - 500 tested frequencies. Next, we determined 
the maximum peak in each power spectrum and show the results 
in Pig. lA.ll In the left panel, the results are shown for the case 
of a mean count rate of 20 counts per second. The light grey 
line indicates the evolution of random maximum peak powers 
with trial number, while the thick black line shows the same 
data sorted by peak powers. To estimate the probability of a 
random encounter at a given confidence level, the sorted curve 
can be interpolated. Eor example, for a detection at 99% we 
determined the power value that corresponds to the sorted curve 
at trial number 9000 of the black line, in this case, yielding 
z = 9.19. We extracted the thresholds for four confidence levels 
and included the results in the right legend of the left panel for 
synthetic light curves with a mean count rate of 20 counts per 
second. We performed the same simulations for mean count 
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Fig. A.2. Top: Comparison of 99% detection thresholds obtained 
with simulations as in Fig. lA. B using different oversampling fac¬ 
tors 16 {Nf - 500), 3 (Nf = 90), and 1 (Nf - 30) with the 
number of tested frequencies Nf given in the legend. The dotted 
lines are the thresholds obtained from Eq. nA.21 i. and they agree 
best with the simulation results without oversampling, while for 
higher oversampling rates, the required power for a 99% detec¬ 
tion is much lower. Bottom: shift in power as a function of de¬ 
tection probability, illustrating that it is not constant. When de¬ 
termining detection probabilities, we correct for oversampling 
by interpolating these simulations. 
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rates ranging from 0.1-350 counts per second and show the 
results in the right panel of Fig. lA.ll The simulation results 
are shown with t he black bullet points, and the predictions 
by IScargld (Il982h . reformulated with Eq. (IA.2b . are included 
for the same confidence levels as the blue dotted lines. No 
systematic trend with count rate can be identified, but all 
simulated results are well below the corresponding theoretically 
predicted thresholds. 

To understand the discrepancy, we repeated the simula¬ 
tions with two additional oversampling factors, thus testing 
90 frequencies (oversampling by 3) and 30 frequencies (no 
oversampling). In Fig. IA.2I we show the results only for the 
99% detection thresholds for these three oversampling factors 
with black, red, and blue as labelled in the upper right panel. 
In our case, using a period range 25-lOOs (df = 0.03 Hz) 
and light curve duration 1000s, the given values Nf need to be 
divided by 30 to derive the oversampling factor. The dotted 
lines are the thresholds calculated from Nf using Eq. (IA.21 i. and 
the dotted and solid lines agree if no oversampling is applied. 
The discrepancies can therefore solely be attributed to the 
oversampling. In the bottom panel of Fig. IA.2I the difference 
between powers with and without oversampling is shown for 
the same three oversampling factors as a function of detection 
probability, from which it can be seen that the discrepancy 
depends on the confidence level. Therefore, for this work, we 
converted observed power values by interpolating look-up tables 
computed from the simulations instead of using Eq. (IA.21 l to 
account for the oversampling effects. 
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An additional effect that needs to be taken into account is 
multiple testing when probing a series of time se gments. We 
applied the conservativ e Bonferroni correction dBonferronil 
Il936tlDunnllT95^Il96lh to the detection thresholds by dividing 
the p value by the number of time segments. For example, 
if we hnd a single peak in a series of 50 power spectra and 
wish to know whether this is signihcant at 99% confidence, we 
divide the p-value of 0.01 by 50 for the test. A 99% detection 
is thus only achieved if the power value exceeds a value that 
corresponds to X = (1 - 0.01 /50) = 99.98%. 


18 







